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Abstract. Multi-epoch multi-wavelength spectroscopic observations of photospheric-phase type II supernovae 
(SN) provide information on massive-star progenitor properties, the core-collapse mechanism, and distances in the 
Universe. Following successes of recent endeavors (Dessart & Hillier 2005ab) with the non-LTE model atmosphere 
code GMFGEN (Hillier & Miller 1998), we present a detailed quantitative spectroscopic analysis of the type II 
SN1999em and, using the Expanding Photosphere Method (EPM) or synthetic fits to observed spectra, a la Baron 
et al. (2004), we estimate its distance. Selecting eight epochs, which cover the first 38 days after discovery, we obtain 
satisfactory fits to optical spectroscopic observations of SN1999em (including the UV and near-IR ranges when 
available). We use the same iron-group metal content for the ejecta, the same power-law density distribution 
(with exponent n = 10 — 12), and a Hubble-velocity law at all times. We adopt a H/He/C/N/0 abundance 
pattern compatible with CNO-cycle equilibrium values for a RSG/BSG progenitor, with C/O enhanced and N 
depleted at later times. The overall evolution of the spectral energy distribution, whose peak shifts to longer 
wavelengths as time progresses, reflects the steady temperature/ionization-level decrease of the ejecta, associated 
non- linearly with a dramatic shift to ions with stronger line-blocking powers in the UV and optical (Fell, Till). 
In the parameter space investigated, CMFGEN is very sensitive and provides photospheric temperatures and 
velocities, reddenings, and the H/He abundance ratio with an accuracy of ± 500 K, ± 10%, 0.05 and 50%, 
respectively. Following Leonard et al. (2002), and their use of correction factors from Hamuy et al. (2001), we 
estimate an EPM distance to SN1999em that also falls 30% short of the Cepheid distance of 11.7 Mpc to its host 
galaxy NGC 1637 (Leonard et al. 2003). However, using the systematically higher correction factors of Dessart 
& Hillier (2005b) removes the discrepancy. A significant scatter, arising primarily from errors in the correction 
factors and derived temperatures, is seen in distances derived using different band passes. However, adopting 
both correction factors and corresponding color-temperatures from tailored models to each observation leads to 
a good agreement between distance estimates obtained from different band passes. The need for detailed model 
computations thus defeats the appeal and simplicity of the original EPM method, which uses tabulated correction 
factors and broadband fluxes, for distance determinations. However, detailed fits to SN optical spectra, based 
on tailored models for individual SN observations, offers a promising approach to obtaining accurate distances, 
either through the EPM or via the technique of Baron et al. (2004). Our best distance-estimate to SN1999em is 
11.5±1.0 Mpc. We note that to achieve 10-20% accuracy in such distance estimates requires multiple observations, 
covering preferentially a range of early epochs preceding the hydrogen-recombination phase. 
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1. Introduction 

Stellar explosions, due to their large luminosity, repre- 
sent an excellent probe of the distant, and thus young. 
Universe. Short Gamma Ray Bursts, identified with the 
death of Wolf-Rayet stars (MacFadyen & Woosley 1999; 
MacFadyen et al. 2001), rank top of the prospective list, 
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but currently, these phenomena are not suitably under- 
stood for distance calibration purposes. Down on the ener- 
getics ladder, supernovae explosions, discovered as bursts 
of optical light, have thus far been a viable route. The 
inherent uniformity of Chandrasekhar-mass white-dwarfs, 
associated with type la SN, has motivated their nomina- 
tion as standard candles. For the accuracy level sought, 
however, dust obscuration, metallicity effects, and (un- 
known) progenitor differences might invalidate such an as- 
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sumption. Although typically ten times less luminous than 
type la, type II SN, resulting from the collapse of the Fe- 
corc of a massive star, represent an alternative. By assiim- 
ing their ejecta expand homologously and radiate a black- 
body spectrum at a well-defined temperature, Kirshner 
& Kwan (1974) designed the Expanding Photosphere 
Method (EPM), a variant of the Baade (1926) method 
for variable stars, to constrain the distance to a SN. The 
agreement between the EPM-distancc to SN1987A and 
that obtained, e.g., with LMC Cepheids, warranted appli- 
cation to a bigger sample, stretching into the Hubble flow 
and thus providing some independent constraint on the 
Hubble constant (Schmidt et al. 1992). But the 30-50% 
lower value of the EPM-distance to SNlQQQem (Leonard 
et al. 2002, L02; Hamuy et al. 2001, HOI) compared to the 
subsequently obtained Cepheid-distance to its host galaxy 
(NGC 1637; Leonard et al. 2003) led to question the reli- 
ability of the EPM. Using the Spectral-fitting Expanding 
Atmosphere Method (SEAM), Baron et al. (2004, B04) 
found a good agreement between their inferred distance 
to SN1999em and the Cepheid-distance to NGC 1637, and 
suggested that correction factors, used to approximate the 
Spectral Energy Distribution (SED) with that of a single- 
temperature blackbody, were too small at later epochs. 

In this paper, we study the various ingredients enter- 
ing the EPM to hunt for the origin of the above discrep- 
ancy, surprising because both methods share a number of 
assumptions or uncertainties: they both assume homolo- 
gous expansion, are both affected by uncertainties in out- 
flow expansion velocity or photometric measurements, and 
rely on multi-epoch observations to constrain the poorly- 
known time of explosion. Besides the generic 20% scatter 
of correction factors, the 20-50% lower values found by 
Eastman et al. (1996, E96) compared to those of Dessart 
& Hillier (2005b, Paper II) could partly or fully explain the 
underestimate of the distance to NGC 1637 and SN1999em 
based on the EPM (HOI, L02). Our approach to determin- 
ing distances to type II SN requires a detailed modeling of 
multi-epoch spectroscopic observations; we thus also infer 
progenitor and SN outflow properties. 

In Dessart & HilHer (2005a, Paper I), we presented 
promising first analyses of photospheric-phase type II SN 
with CMFGEN (HiUier & Miller 1998), a model atmo- 
sphere code so far devoted to stiidics of massive star winds 
(Najarro et al. 1997, Hillier & Miller 1999, Dessart et al. 
2000, Crowther et al. 2002, Martins et al. 2004). To step 
progressively through the various challenges posed by SN 
spectroscopic modeling, we first concentrate on type II 
SN, known to exhibit modest, and at early times negligi- 
ble, amounts of nucleosynthesised metals during the ex- 
plosion, thereby allowing us to neglect the deposition of 
energy from radioactive decay. Together with a Hubble ve- 
locity law, we implemented a simple power law density dis- 
tribution p{r) = po{Ro/r)'^, where po and Rq are the base 
density and radius, and n is the density exponent. Using 
spectroscopic observations of SN1987A and SN1999em, se- 
lected at times when hydrogen is partially or fully ion- 
ized, we found that both objects could be well modeled 



by adopting standard CNO-cycle equilibrium abundances 
typical of a supergiant massive star progenitor. Two op- 
tical features in early-time spectra, previously associated 
with peculiar line emission/ absorption processes and/or 
outflow inhomogeneities lying high above the photosphere 
(Baron et al. 2000) , were unambiguously identified as N II- 
multiplet lines. The strength of the features suggested a ni- 
trogen enrichment, at the corresponding epoch, of a factor 
of a few compared to cosmic for SN1999em, but somewhat 
less for SN1987A. Hei5875A could also be weh repro- 
duced with a modest helium enrichment, in contrast with 
former difficulties associated with that species (Eastman 
& Kirshner 1989). Additionally, optical line profiles, and 
in particular their strong and persistent blueshift, were 
found to place a strong constraint on the density-gradient 
in the vicinity of the photosphere, requiring n ~ 10. For 
these test cases, CMFGEN fitted the UV, the optical and 
the near-IR (when available) with an accuracy generally 
better than a few tens of percent. 

Overall, this high-level of agreement between our syn- 
thetic SEDs and spectroscopic observations of type II SN 
during their photospheric phase validates such a detailed 
approach with CMFGEN. It allows us to draw conclu- 
sions on the progenitor properties and perhaps more im- 
portantly on distances to the host galaxy. 

In Dessart & Hillier (2005b), we reviewed the various 
items entering the EPM, to complement the only other 
study of the kind done by Eastman et al. (1996, E96). 
CMFGEN yields correction factors that are systematically 
^20% larger than those of E96, but confirms the general 
properties with effective- or blackbody color-temperature, 
and photospheric density. We also stressed that correction 
factors do not only correct for dilution by electron scat- 
tering of the SED set at the thermalization depth, but, 
within the standard domain of application of the EPM, 
also correct for the corrupting effects of line emission and 
absorption on the continuum flux within a given band- 
pass. Following the appearance of metal liric^ blanketing 
below Tphot ^ 8000K, such line effects dominate in the 
S-band (attributable to Fe ii and Ti ii) and in the /-band 
(Ca ii), making the SED incompatible with that of a black- 
body. Due to the steep increase of correction factors at low 
color-temperatures, modulations of line-strengths (having 
nothing to do with ionization but resulting from differ- 
ent metal abundances or expansion velocity) can affect 
inferred color-temperatures and generate substantial er- 
rors. Moreover, associated optical-depth enhancements in 
the corresponding spectral regions make the photosphere 
more extended than in the nearby line-free regions, inval- 
idating the uniqueness of the photosphere radius required 
by the Baade method. These conclusions thus argue for 
the use of the EPM at early-times, when the outflow is 
fully ionized and metal-line blanketing conflned to the UV 
region. Finally, although the velocity at maximum absorp- 
tion Umax in optical P-Cygni profiles corresponds to the 
photospheric velocity tiphot to within few tens of percents, 
we showed that Wmax can over- as well as under-estimate 
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phot 5 i-e. Vmax does not necessarily exceed Wphot (in abso- 
lute terms). 

In this paper, we analyze a high-quality photometric 
and spectroscopic dataset for SN1999em with CMFGEN 
and infer its distance, first with the EPM in combina- 
tion with various sources for correction factors, and then 
using CMFGEN fits to each epoch (an approach similar 
in spirit to the SEAM of B04). In Sect. [3 we present 
our selection of spectro-photometric observations, accord- 
ing to the criteria outlined above. In Sect. 01 we show 
the results of the quantitative spectroscopic modeling for 
each of the selected dates, drawing information on the 
outflow kinematics, chemistry, ionization. In Sect. ^ we 
compute the distance to SN1999em. In Sect. 14.11 using 
the standard minimization procedure and prescriptions 
for correction factors from HOI, the EPM-distance shows 
a large scatter between selected band passes, is compa- 
rable to that found by L02 and HOI and thus always 
lower than the Cepheid-distance to NGC 1637. When we 
use the correction factors from Paper II, distance esti- 
mates are in better agreement with the Cepheid-distance, 
but there is still a relatively large scatter between dif- 
ferent filter sets. A similar scatter occurs when we use 
blackbody color-temperatures extracted from the detailed 
models with the correction factors tabulated in Paper II. 
As expected, consistent distance estimates are obtained 
when we use the model blackbody color-temperatures to- 
gether with the model correction factors. For comparison, 
we compute the SEAM-distance to SN1999em in Sect.lO 
and corroborate the EPM-distance determination based 
on CMFGEN models or that of B04. In Sect.El we draw 
our conclusions and discuss future prospects. 

2. Spectrophotometric datasets of and reddening 
to SN1999em 

No neutrino detection was associated with the SN1999cm 
event and thus no precise dating of the core-collapse of the 
progenitor exists. Inferences, based on the light curve ap- 
pearance or the EPM/SEAM, suggest the explosion took 
place around the 25th of October 1999, with an uncer- 
tainty of ± 2 days (HOI, L02, B04). It was absent in im- 
ages of the host galaxy NGC 1637 taken on October 20.45 
UT, while first detection of the SN was made on October 
29th (JD 2451480.94; Li 1999). The resemblance of the 
spectrum taken on October 30th with the spectrum of 
SN1987A taken on the 24th of February 1987, indicates 
that SN1999em was indeed caught at a very early phase 
of evolution. 

In this work, we use two main sources of spectroscopic 
observations to model the evolution of SN1999em, from 
the first spectrum taken on the 30th of October 1999 un- 
til the last one taken on the 5th of December 1999, by 
which time metal line-blanketing in the optical spectral 
range is strong. From HOI, we use optical spectra taken 
on the 30th of October, and the 3rd, 9th, 14th, and 19th 
of November. In addition, we use the near-IR observa- 
tions of HOI on the 2nd and 19th of November. From the 



SUSPECT^ archive, we extract optical data of L02 for 
the 1st and 5th of November, and the 5th of December. 
Finally, we use HST-STIS observations taken on the 5th 
of November, kindly provided by Baron (priv. comm.), 
and covering the wavelength range from 1140A to 5710A, 
overlapping in the red with observations of L02. On that 
date, we merge the two datasets to have a full coverage of 
the UV and optical out to 7500A. We exclude from our 
sample the observations of the 22nd and 24th of November 
(L02) because of dubious absolute flux level, curious rel- 
ative flux distribution at and beyond Ha (compared to 
spectra taken at embracing dates; see discussion in L02). 
Hence, our sample contains eight observations for which 
photometric measurements have been obtained by L02, co- 
incident to within a few hours to spectroscopic ones. On 
day 38, we adopt the magnitudes of L02 and combine these 
measurements with the spectrum taken on the previous 
day, i.e., the 5th of December, justified by the very slow 
spectroscopic evolution at such late times. In Tabled we 
provide a log of the observations used in this work. Despite 
HOl's photometric measurements over a large time-span 
in the Z , J, H , and K bands, we limit our investigation 
to the optical band passes B, V, and /, and sets {B, V}, 
{B,V,I}, and {V, 1} (note one exception with the Z-band 
in Sect. 14.21 see HOI for a presentation of this Z-band). 

Our modeling procedure blueshifts the spectroscopic 
observations of SN1999em by 770 km s"^ (L02), scales the 
synthetic flux to the Cepheid-distance of 11.7Mpc (L03), 
and finally reddens such a scaled synthetic spectrum using 
the CardeUi et al. (1998) law, adopting Av/E{B -V)^ 
3.1. In Paper I, we illustrated the mingled effects on the 
SED associated with metal line-blanketing, cooling of the 
outflow and reddening. Here, the fine time-sampling and 
the HST-STIS data can only be fitted with CMFGEN by 
adopting a low reddening value of 0.1±0.05, in agreement 
with previous studies (L02, Baron et al. 2000, B04; see 
also Sect. 13.3^ . Note that the UV data is indeed a key as- 
set for reddening determinations since the SED at longer 
wavelength becomes less and less sensitive to reddening 
modulations. In the context of the EPM, the extinction 
uncertainty translates directly into an uncertainty in the 
distance modulus. Fortunately, due to the low ISM column 
density to SN1999em, this introduces modest absolute er- 
rors on the distance (see Sect. 14. 2|) . 

To highlight line-emission contributions, we always 
plot the synthetic spectrum obtained by including all 
opacity sources (red curve) and that accounting only for 
continuum processes (blue curve). Finally, the adopted 
model luminosities for each date match the absolute ob- 
served flux level (for the assumed Cepheid-distance) to 
within a factor of approximately 2-3 (our initial adoption 
of the lower distance estimate by HOI /L02 when modelling 
SN1999em introduced a systematic underestimate of the 
luminosities by a factor of ^(11.7/8)^ ^^-"2.1); the synthetic 
SED shows only a modest sensitivity to luminosity scal- 
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Table 1. Log of the observations used and described 
in Sect. HJ Photometric measurements, with correspond- 
ing dates, are reproduced from L02 (columns 2-5). 
Spectroscopic datasets are referenced according to spec- 
tral range, i.e. UV, optical (Opt.) and near-IR (nIR). Note 
that in Sect. ^ the photometric data on day 38 is com- 
bined with the spectrocopic data on day 37. 



Date 


Day" 


Photometry 


b 


Spectroscopy 


(1999) 


B 


V 


I 


UV Opt. nIR 


30 Oct 


1.0 


13.87 


13.87 


13.65 


1 


01 Nov 


3.0 


13.80 


13.79 


13.56 


2 


02 Nov 


4.0 








1 


03 Nov 


5.0 


13.85 


13.79 


13.54 


1 


05 Nov 


7.0 


13.92 


13.84 


13.53 


3 2,3 


09 Nov 


11.0 


14.02 


13.84 


13.48 


1 


14 Nov 


16.0 


14.25 


13.81 


13.44 


1 


19 Nov 


21.0 


14.47 


13.86 


13.40 


1 1 


06 Dec 


38.0 


14.94 


13.93 


13.29 


2 



a: days since JD 2,451,480.94; b: measurements are taken 
from 1: HOI; 2: L02; 3: Baron et al. (2000). 



ings of that magnitude, provided the outflow ionization is 
preserved (Paper I, E96). 

In Sects. 14.114.21 we convert synthetic and blackbody 
fluxes into magnitudes using the filter transmission func- 
tions and zero-points of HOI and Hamuy (priv. comm.), 
adopting either the value at each filter's effective wave- 
length {Ab = A.05E{B-V) and Aj = 1ME{B-V)), or 
convolving the synthetic SED with the filter transmission 
function. Resulting magnitudes from the two approaches 
differ by no more than 0.005 mag, thus much less than 
other uncertainties involved in the EPM or SEAM. 

3. Quantitative spectroscopic analysis of 
SN1999em 

In this section, we perform a quantitative spectroscopic 
analysis of SN1999em, at eight dates sampling its early 
photospheric-phase evolution (see above). In Sect. 13. ll we 
review the standard model parameters and assumptions. 
We then detail our results for each observation in Sect. 13.21 
Finally, in Sect. 13.31 we discuss typical uncertainties asso- 
ciated with each of the inferred model parameters. 

3.1. Model Presentation 

The standard set of parameters needed by CMFGEN were 
discussed thoroughly in Paper I. A key asset of CMFGEN 
is the explicit treatment of line-blanketing, justifying a de- 
tailed description of the atomic structure for a large num- 
ber of species. Models discussed below include Hi (30,20), 
Hei (51,40), Ci (63,33), Cii (59,32), Cm (20,12), Civ 
(14,9), Ni (104,44), Nil (41,23), Niii (8,8), Oi (75,23), 
On (111,30), Oiii (46,26), Neii (242,42), Nai (71,22), 
Nan (35,21), Mgii (65,22), Alii (44,26), Aim (45,17), 



Sin (59,31), Sim (51,27), Sn (324,56), Sin (98,48), Can 
(77,21), Tin (152,37), Tim (206,33), Cm (196,28), Crin 
(145,30), Mnn (97,25), Mnni (175,30), Fei (136,44), Fen 
(309,116), Feni (477,61), Feiv (282,50), Con (144,34), 
Com (283,41), Nin (93,19), Nini (67,15), each parenthe- 
sis containing the number of full- and super-levels (see 
Hillier & Miller 1998 for details). At epochs when a given 
species offers a negligible opacity to the radiation field at 
all depths in the outflow, we exclude it from the atomic 
dataset, thereby reducing the memory requirements and 
accelerating the convergence of CMFGEN. Under typical 
type II SN conditions and when all the above species are 
included, CMFGEN uses up to 3.5GB of RAM and may 
take a week to reduce all level population fluctuations, 
from one iteration to the next, to less than 1%. 

As discussed below, the main component controlling 
the evolution of type II SN spectra during their photo- 
spheric phase is the steady decrease of the outflow tem- 
perature, in parallel with the reduction of the outflow 
ionization. As time proceeds, opacity sources change, and 
the opacity may increase (e.g., when Fen dominates over 
Fe in). In contrast to the large changes in effective temper- 
ature, we flnd only small and specific changes in chemical 
composition as time progresses through the first 40 days 
past discovery. No metal-abundance variations associated 
with explosive nucleosynthesis are found, so that all spec- 
tra are fitted with a unique metallicity, adopted to be solar 
(L02). For hydrogen and helium, we reproduce all epochs 
with a unique value, H/He = 5 (by number), although our 
inability to fit the Balmer lines at late times makes our re- 
sult reliable only up to the hydrogen-recombination phase 
(first twenty days) . The nitrogen abundance is constrained 
from the Nil optical lines (Paper I) and requires over- 
solar abundance; we adopt N/He = 6.8 xlO"^ (Prantzos 
et al. 1986). C and O abundances are diflicult to con- 
strain: at early times, since no lines from these species 
are observed, we adopt C and O abundances compatible 
with the N/He above and CNO-cycle equilibrium values, 
with C/He = 1.7 xlO"'' and 0/He = 10"'' by number. At 
later times, however, C i and O i lines are unambiguously 
observed in the 7000-12 000 A region, requiring enhance- 
ments compared to CNO equilibrium values, of a factor of 
10 for carbon and 100 for oxygen fSect. lTTH|l . Finally, the 
Na 1 5890A doublet can only be fitted with a factor of four 
over-abundance compared to solar, in agreement with the 
CNO-cycle equilibrium values of Prantzos et al. (1986), 
i.e. A^Na = 1-38 X 10^**. Overall, we note that the con- 
spicuous spectroscopic evolution during the photospheric 
phase is not a refiection of outflow changes in chemical 
composition but rather in ionization/temperature. 

The density distribution of the outflow is constrained 
via the density exponent n. In Papers I and II, we il- 
lustrated, from different perspectives, the spatial conflne- 
ment of line and continuum formation regions to a few 
density scale heights about the photosphere rather than 
over the entire outflow - the Can set of lines at around 
8500A is a noticeable exception (e.g.. Fig. 13., Paper I). 
The resulting line-profile shapes require values n ^ 10 to 
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Table 2. Summary of model results presented in Sect. 13.21 Time reference is in days since JD 2,451,480.94. On the 
right hand side, we also give the key EPM quantities for each model and date, i.e., blackbody color temperature and 
corresponding correction factor for each bandpass, {B, V}, {B, V, I}, and {V, 1} (Sect. 01 Paper II). 



rim/ 
uay 


r 


-t phot 


Jtphot 




Pphot 


Tl 


-t BV 


S,BV 


J BVI 


£,BVI 


-t VI 






I'm'* r 




liu cm I 


I Kill S ) 


l^iu gem ) 
















1.0 


15.0 


13.38 


6.2 


10,430 


4.3 


12 


16.02 


0.462 


14.96 


0.50 


14.39 


0.52 


3.0 


15.0 


12.59 


6.8 


8,940 


3.3 


10 


15.96 


0.430 


15.03 


0.46 


14.23 


0.49 


4.0 


15.0 


12.59 


6.8 


8,940 


3.3 


10 


15.96 


0.430 


15.03 


0.46 


14.23 


0.49 


5.0 


9.0 


10.63 


6.7 


8,840 


3.6 


10 


14.23 


0.443 


13.20 


0.49 


12.53 


0.51 


7.0 


5.0 


9.20 


6.64 


8,750 


4.1 


10 


12.95 


0.443 


11.95 


0.49 


11.28 


0.53 


11.0 


2.5 


8.04 


6.65 


7,960 


4.0 


10 


10.70 


0.469 


10.25 


0.50 


9.83 


0.53 


16.0 


1.5 


6.80 


6.15 


6,350 


8.7 


10 


8.17 


0.646 


8.84 


0.56 


9.55 


0.51 


21.0 


3.75 


6.26 


12.90 


5,530 


7.4 


10 


7.11 


0.685 


7.85 


0.56 


8.74 


0.48 


38.0 


2.0 


5.92 


10.90 


3,400 


38.8 


10 


5.54 


1.072 


6.34 


0.75 


7.24 


0.61 



reproduce the blueshift as well as both the modest absorp- 
tion and emission observed in Balmer lines. Comparable 
values are obtained in other spectroscopic studies, e.g. 
Eastman & Kirshner (1989), or independent radiation- 
hydrodynamical simulations of core-collapse (Arnett 1988; 
Woosley 1988; Ensman & Burrows 1992). 

We quote numerous model properties at the photo- 
sphere, which corresponds here to the location where the 
continuum optical depth, integrated inwards from the 
outer grid radius, is 2/3. We constrain the expansion, or, 
equivalently, the photospheric velocity by adjusting the 
velocity at the inner boundary; we then adopt the value 
that permits a satisfactory fit to profile shapes for a wide 
range of optical line diagnostics, a procedure we favor over 
the direct and more precise, but less meaningful, measure- 
ments on individual observed line profiles. Further discus- 
sion on this is provided in Sect. 5 of Paper II, as well as 
below, in Sect. 13.31 

3.2. Tailored Analysis with CMFGEN 

In this section we model the data for our selected obser- 
vation dates (Sect. EJ, provide synthetic fits to observed 
SEDs (Figs.dHl), tabulate the derived model param- 
eters (Table As described in the preceding section, the 
outfiow temperature, modulated via the base luminosity 
or radius, represents the essential change between con- 
secutive dates. The electron temperature at the photo- 
sphere, Tphot, accounts well for this change, and thus con- 
stitutes the center of the discussion below. To infer the 
photospheric velocity, many models were calculated: we 
adopted the photospheric velocity corresponding to that 
model which provided the best fit to the majority of ob- 
served line profiles. A detailed discussion on the uncer- 
tainties in the derived photospheric velocity is provided in 
Sect.ESSl 

3.2.1. Observation of the 30th of October 1999 

In the top panel of Fig. ^ we present a fit of the red- 
dened scaled synthetic spectrum (full: thick red; contin- 



SN1999em - 30/10/99 - E(B-V) = 0.1 











— Observed Spectrum 






— Total Reddened Synthetic Spectrum 






— Continuum Reddened Synthetic Spectrum 





4000 5000 6000 7000 8000 9000 

A (A) 

SN1999em - 30/10/99 - E(B-V) = 0,1 




4000 5000 0000 7000 8000 9000 

A (A) 



Fig. 1. Top: fit of the full (red) and continuum (blue) 
synthetic spectra to observations (black) of SN1999em 
taken on the 30th of October 1999 (Day 1; HOI). Model 
parameters are: = 1.5 x lO^L©, Tphot = 13 380K, 
-Rphot = 6.2 X 10^'^cm, Uphot =10430kms~\ pphot = 
4.3 X 10-i'*gcm"3^ g^j^^ ^ ^ 12 The synthetic fiux, 17% 
lower than observed, is re-normalized at 8200A. Bottom: 
same as above, but this time overplotting reddened syn- 
thetic spectra computed by treating bound-bound transi- 
tions of a single species: Hi, He i, Nil, and Feiii. Contrary 
to Baron et al. (2000), we attribute the two features just 
blue-ward of H/3 and He 1 5875A to N ii lines. (This figure 
is available in color in the electronic version.) 



uum: thin blue) to the observations taken on the 30th 
of October 1999 (black; HOI). The model parameters are: 
= 1.5 xIO^Lq, Tphot = 13380K, i?phot = 6.2xl0i''cm, 
Wphot = 10 430kms~\ pphot = 4.3 x lO^^'^gcm"^, and 
n = 12. The flux, underestimated by 17%, is renormalised 
at 8200A. The overall fit quality is good, although the fiux 
around 3800A is underestimated. Including more species 
or bigger model atoms makes no difference; increasing 
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Iphot helps but then provides a poor match to the ob- 
served flux elsewhere, e.g., below 3500A. A poor flux cali- 
bration in this spectral region could be at the origin of this 
discrepancy. The fit quality to Balmer lines is satisfactory, 
although the extent of profile troughs is underestimated. 
At such early times, our neglect of relativistic terms in the 
radiative transfer equation may contribute to such discrep- 
ancies (Hauschildt et al. 1991; Jeffery 1993). Helium lines, 
and most notably He 1 5875A, are well reproduced. In the 
bottom panel of Fig. ^ we show synthetic spectrum fits 
(thin colored-lines) to observations, but treating chemical 
species in isolation. As emphasized in Paper I, we predict 
N II (orange line) features in the blue wing of Hei5875A, 
resulting from the 3d-3p multiplet around 5470A and the 
3p-3s multiplet around 5670A, as well as in the blue wing 
of H/3, resulting from other Nil multiplets around 4600A. 
Finally, notice how the continuum lies in general well be- 
low the full synthetic flux distribution level, even in spec- 
tral regions that look essentially featureless, revealing the 
ubiquitous presence of line emission (see, e.g., the subtle 
contribution from Fe ill - colored in red - to the red wing 
of both H/3 and Hei5875A, or that of Hei6678A in the 
red wing of Ha). 

3.2.2. Observations of the lst-2nd of November 1999 

In Fig. [3 we reproduce Fig. for the observations of the 
1st (optical, top panel; L02) and the 2nd (near-IR, bot- 
tom panel; HOI) of November 1999. The fuU (thick red) 
and continuum (thin blue) synthetic flux distributions un- 
derestimate by 30% the observed flux and are thus re- 
normalized, at 6900A, to fit both spectral ranges (with 
an additional 8% flux adjustment between optical and 
near-IR). The model parameters, identical for both spec- 
tral ranges, are: = 1.5 x IO^Lq, Tphot = 12 590K, 
-Rphot = 6.8 X lO^'^cm, Uphot = 8940kms"\ Pphot = 
3.3 X 10-"gcm-3, and n = 10. 

In the optical, there is little change with observations 
taken on the 30th of October, reflected in the similar- 
ity of model parameters for both dates. We obtain good 
fits to Balmer as well as Paschen lines, i.e., P7 1.094 /im, 
P/31.28/im, and Pa 1.875 /xm. Note, however, that profile 
troughs are sometimes predicted too narrow (e.g.. Ha and 
H/3). He I lines are also well reproduced, both in the opti- 
cal at 5875 A and in the near-IR at 1.083 /im (which over- 
laps with P7 but contributes most of the flux seen at this 
date in the 1.08 /im feature). Note the broad Nil feature 
blueward of Hei5875A, although its absorption strength 
is somewhat overestimated (see Sect. 13. 3| . 

3.2.3. Observation of the 3rd of November 1999 

In Fig. 13 we present a flt of the full (red) and contin- 
uum (blue) synthetic SED to the observed SED of the 
3rd of November 1999 (black; HOI). The model has the 
following properties: = 9 x lO**^©, Tphot = 10 630K, 
-Rphot = 6.7 X lO^'^cm, Uphot = 8840kms"\ Pphot = 
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Fig. 2. Fit of the full (red) and continuum (blue) syn- 
thetic spectra to observations (black) of SN1999em taken 
on the 1st (optical, L02, day 3; top) and the 2nd (near- 
IR, HOI, day 4; bottom) of November 1999. In the bot- 
tom panel, we include the synthetic spectra obtained 
by restricting bound-bound transitions to those of H i 
(blue) or Hei (green). The model parameters are: — 

1.5 X lO^L©, Tphot = 12 590K, i?phot = 6.8 x lO^^cm, 
Vphot — 8940 km s^^, Pphot — 3.3 x 10^^'^gcm^'^ and 
n = 10. The model flux, which underestimates the ob- 
servations by 30%, is renormalized to the observed value 
at 6900A. (This flgure is available in color in the electronic 
version.) 

3.6 X 10"^"* gcm~'^, and n = 10. The model flux, under- 
estimated by 50%, is renormalised to the observed value 
at 5400A. The fit quality is good. Increasing the outflow 
temperature (or ionization) results in a better flt of the 
overall shape of the spectrum but leads to an overestimate 
of the strength of He 1 5875A and that of the N ii 5400A 
feature. The bump in the red wing of H/3, due primar- 
ily to Fe III lines, is underestimated, suggesting perhaps 
a slightly higher environmental metallicity (Paper I). As 
noted in the preamble to the section, we use, by default, a 
solar mixture of metals. Finally, note how reduced the con- 
tinuum flux level is compared to the full SED, highlighting 
the ubiquitous presence of broad and weak "background" 
lines. 

3.2.4. Observation of the 5th of November 1999 

In Fig. 01 we show a full (red) and continuum (blue) syn- 
thetic flt to observations taken on the 5th of November 
1999 (L02), the optical range being complemented with 
UV HST data (Baron et al. 2000; see Sect.©. The model 
has the following parameters: T* = 5 x lO^L©, Tphot — 
9200 K, i?phot = 6.64 x lO^^^cm, Wphot = 8750 km s'^, 
Pphot — 4.1 X 10~^'*gcm~'^, and n = 10. The model flux, 
which underestimates the observations by a factor of two, 
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confirmation, of the reddening obtained usually at other 
(relatively late) times from optical data. 



Fig. 3. Synthetic fits (color) to observations (black) taken 
on November 3rd (HOI; day 5). The model has the fol- 
lowing parameters: = 9 x IO^Lq, Tphot = 10 630K, 
-Rphot = 6.7 X 10^** cm, Uphot = 8840kms~\ Pphot = 
3.6 X 10^^"' gcm"-^, and n = 10. The model flux, which 
underestimates the observations by 50%, is renormalized 
to the observed value at 5400A. (This figure is available 
in color in the electronic version.) 



is renormalized to the observed value at 7500A. Both ob- 
servation and model were discussed in detail in Paper I 
(Sect. 3.2), emphasizing the complex and primary role 
played by metal line-blanketing in the UV range (see lower 
panel of Fig. 3, Paper I). 

The fit quality is good but there are noticeable dis- 
crepancies. The optical flux is somewhat underestimated, 
while the UV flux is overestimated in a number of wave- 
length regions. We overestimate the flux in the 2800A re- 
gion, which, additionally to Fe ii, is strongly influenced 
by the multiplet of Mgii, composed of Mgii (3p-3s) and 
Mgii (4s-3p). Increasing/decreasing the Magnesium abun- 
dance by a factor of two introduces flux changes only 
at the 10% level, thus insufficient to resolve the discrep- 
ancy. In Paper I, we showed how a global increase of 
the pristine metal abundance leads to a better match, 
e.g., of the strength of the iron feature at '--^SOOOA, giv- 
ing some support for a solar or over-solar metallicity of 
the SN1999em environment; simultaneously, the increased 
metal-line blanketing in the UV reduces the emergent UV 
flux to more suitable levels. Alternatively, lowering the 
effective temperature increases line-blocking through the 
shift from Feiii to Fell opacity sources, also matching 
the observed UV flux better; in this situation, however, 
the slope of the optical synthetic SED is then too flat. 
We discuss further these issues in Sect. 13.31 but already 
here, we see how coupled and subtle the effects of line- 
blanketing and outflow temperature/ionization are. Note 
that the photospheric temperature of 9200 K found here 
is in good agreement with the value of 9000 K proposed 
for that date by Baron et al. (2000). 

Overall, having access to UV data during the early 
evolution is an asset, as it allows the determination, or the 
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Fig. 4. Full (red) and continuum (blue) synthetic flts to 
observations (black) taken on November 5th (UV range: 
Baron et al. 2000; optical range: L02, day 7). The model 
has the following parameters: L, = 5 x lO^i©, Tphot = 
9200K, i?phot = 6.64 x 10"cm, Wphot = 8750 kms^^, 
Ppjjot = 4.1 X lO^^^gcm"'^, and ?i = 10. The model flux, 
which underestimates the observations by a factor of two, 
is renormalized to the observed value at 7500A. The dip in 
the continuum spectrum near 2200 A is due to the effects 
of interstellar extinction. (This flgure is available in color 
in the electronic version.) 



3.2.5. Observation of the 9th of November 1999 

We show in Fig. a synthetic fit to the observations of 
SN1999em on the 9th of November 1999 (HOI; day 11). 
The model parameters are the following: i» = 2.5 x 
IO^Lq, Tphot = 8040K, ii-phot = 6.65 x lO^^cm, Wphot = 
7960kms-\ pphot = 4 x IQ-i^'gcm-^, and n = 10. The 
model fiux, which underestimates the observations by a 
factor of 2.65, is renormalized to the observed value at 
6OOOA. Starting around this date, metal line blanketing in 
the optical is noticeable, increasing steadily as we progress 
through the photospheric phase. Mingled with the Balmer 
lines (H7, H/3, and Ha), one sees the presence of Can 
at 3800A (as for the previous date). Fell and Tin lines 
at ~4400A, and Fe 11 fines redward of H/3 up to 55OOA. 
He 1 5875A is still observable as a very weak feature, soon 
to be replaced at the same location by the Na 1 5890A dou- 
blet. Fell and Sin are also present around 6300A, creating 
a kink in the Ha absorption trough; it might thus not be 
necessary to invoke a peculiar formation of the Ha line to 
explain this feature (L02). 



3.2.6. Observation of the 14th of November 1999 

In Fig. El we present full (red) and continuum (blue) syn- 
thetic spectrum fits to observations taken on the 14th 



8 



Luc Dessart and D. John Hillier: Type II supernovae spectroscopic modeling 



SN1999em - 09/11/99 - E(B-V) = 0.1 



2,C — 




^DGO 5000 6000 7000 SOOO 9000 

X (A) 



Fig. 5. Synthetic fits (red) to observations (black) taken 
on November 9th (HOI; day 11). The model has the fol- 
lowing parameters: L» — 2.5 x lO^L©, Tphot = 8040 K, 
^phot = 6.65 X 10^"* cm, Uphot = 7960kms~\ pphot = 
4 X 10^"'^"' gcm^'"*, and n — 10. The model flux, which un- 
derestimates the observations by a factor of 2.65, is renor- 
malized to the observed value at 6OOOA. The feature in the 
continuum spectrum near 3600 A arises from the inclusion, 
via level dissolution, of Balmer series bound-bound transi- 
tions near the series limit. (This figure is available in color 
in the electronic version.) 

of November 1999 (HOI; day 16). The model parameters 
are the following: = 1.5 x IO^Lq, Tphot = 6800 K, 
^phot = 6.15 X 10^"* cm, Uphot = 6350kms"\ Pphot = 
8.7x 10~^'*gcm~^, and n = 10. The model flux, which un- 
derestimates the observations by a factor of 3.7, is renor- 
malized to the observed value at 6OOOA. Note that a de- 
tailed discussion of this model has been made in Sect. 3.3 
of Paper L 

A conspicuous feature now appears in the / band, 
attributable to Can 8498-8542-8662 A, consistent with 
the identification of Can 3933-3968A at earlier epochs. 
Hei5875A has also vanished, replaced by a weak feature 
due to Na 1 5890A, well fitted by invoking an enhancement 
of four over the corresponding solar sodium abundance 
(Paper I; Prantzos et al. 1986). The feature at 7700A now 
seems to extend too far to the red to be attributed exclu- 
sively to terrestrial atmospheric absorption, as supposed 
so far, and thus suggests the presence of Oi7700A. In 
the top right corner of Fig. IHl we add, over the relevant 
region, the synthetic SED (green) obtained by enhancing 
the oxygen abundance by a factor of 100 over the standard 
value 0/Hc = 10~*, as required to reproduce the red-side 
of the observed feature. Observations at later epochs, re- 
vealing spectra with features both narrower and better 
centered on the rest-wavelength (see Paper I), will confirm 
this identification and the need for such a significant en- 
hancement of the oxygen abundance over the CNO-cycle 
equilibrium value, characteristic of the deeper layers of the 
progenitor star (Hirshi et al. 2004). Metal line features in 
the optical have also strengthened compared to those of 



the previous date, following the cooling of the outfiow, 
and the recombination of hydrogen and Feiii (to Fell). 
The Ha line profile is not too well fitted: its red-wing ex- 
tent is underestimated as is the P-Cygni trough, which 
is now very pronounced. We instead predict a narrower 
profile shape, with considerable fiUing-in of the trough by 
Fe II and Si 11 lines. As we will see below, this problem only 
gets worse, associated with the infamous underestimate of 
the Ha strength. 
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Fig. 6. Full (red) and continuum (blue) synthetic spectral 
fits to observations (black) taken on November 14th (HOI; 
day 16). The model has the following parameters: — 
1.5 X IO^Lq, Tphot = 6800K, i?phot = 6.15 x lO^^cm, 
Wphot = 6350 kms^^, Pphot = 8.7 x 10^^'^gcm^'^ and n = 
10. The model fiux, which underestimates the observations 
by a factor of 3.7, is renormalized to the observed value 
at 6OOOA. The species identification of line contributions 
is shown in the lower panel of Fig. 5, Paper I. (This figure 
is available in color in the electronic version.) 



3.2.7. Observation of the 19th of November 1999 

We show in Fig.dsynthetic fits (color) to the observations 
of the 19th of November 1999 (HOI; day 21). The model 
fiux is normalized to the observed one at 6IOOA with 
an additional decrease of 15% in the near-IR. This fiux 
shift is likely to result from the time difference between 
the two observations rather than an inadequate slope be- 
tween the optical and the near-IR ranges in the synthetic 
spectrum. The model properties are the following: T, = 
3.75 X IO^Lq, Tphot = 6260 K, i?phot = 12.9 x 10^^ cm, 
Wphot = 5530 kms~^, Pphot — 7.4 x 10^^"* gcm^"^, and 
n = 10. The model flux, which underestimates the ob- 
servations by 10%, is renormalised to the observed value 
at 6IOOA. 

As described in Paper I, we have found that using a 
small turbulent velocity helps maintain signiflcant emis- 
sion in Ha (see next section) . In general, this entails using 
a turbulent velocity of 20kms~^, rather than lOOkms"^ 
used for models where hydrogen is fully or mostly ionized. 
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Note however that this change only delays the onset of the 
discrepancy for the Ha fit, so is introduced here mostly for 
aesthetic reasons. We find that lowering the density expo- 
nent down to 8 and even 6 leads to an increase in the Ha 
strength but also produces a very severe mismatch of the 
Ca II lines, which become both stronger and broader than 
observed. We find that indeed, maintaining a high density 
exponent of 10 leads to better fits to observations, at least 
for the epochs covered here and our model assumptions 
(Paper I). Note also the presence of Fell and Sin lines 
which fill in the Ha trough. 

Omitting this difficulty of reproducing Ha and the 
Paschen series the quality of the fit is again satisfactory, 
with a good match to the strength and shapes of line fea- 
tures due to Fell, Tin, Can. The overall flux distribution 
is very well matched, with a slight overestimate of the flux 
in the blue. Another discrepancy is the slight overestimate 
of the Caii/Fen/Tiii blend around 3800 A. 

We observe a weak feature at ^1.07/im which we at- 
tribute to the C I (3p-3s) multiplet. However, to reproduce 
this weak feature, the carbon abundance had to be in- 
creased by a factor of 10 to C/He = 0.0017, equivalent to 
a mass fraction of 0.0023. Similarly, a number of Ci lines 
over the range 9060"9110A are seen, although a firm iden- 
tification is established with the observation of the 5th of 
December 1999 fSect. 1301^ . Such a carbon abundance is 
compatible with chemical mixing prior to core-collapse, 
inducing a more moderate depletion of C/0 compared 
to N. Accordingly, and following from the previous sec- 
tion, we identify the Oi7700A line (adopting 0/He — 
10~^), as well as another oxygen feature in the near-IR 
at 1.13/xm, standing in the red wing of P7. Note that the 
latter also overlaps with a multiplets of Mg 11 (4p-3d) at 
1. 091-1. 095/im (which represents about 30% of the 1.09yLtm 
feature). 

An observed feature at ~1.03/zm is not predicted by 
CMFGEN, despite the numerous species included; this 
could come from singly-ionized species, or perhaps Sc or 
Ba, not accounted for at present. 

3.2.8. Observation of the 5th of December 1999 

We present synthetic fits (color) to the observations 
(black) taken on the 5th of December 1999 (L02) in Fig.|Hl 
The corresponding model has the following parameters: 
= 2 X lO^i©, Tphot = 5920K, Rph^t = 10.9 x lO" cm, 
Wphot = 3400kms"\ pphot = 38.8 x lO'^'^gcm"^, and 
n — 10. The model flux, which underestimates the obser- 
vations by 10%, is renormalized to the observed value at 
7OOOA. Also, to ease the identification of observed spectral 
features, made more numerous due to the reduced out- 
flow expansion velocity compared to previous dates, we 
show, in the lower panel, rectified spectra obtained when 
only bound-bound transitions of a given ionization species 
(see label) are included, ordered upwards with increasing 
atomic mass. 
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Fig. 7. Synthetic fit (red) to observations (black) taken 
on November 19th (HOI; optical range, day 21) and 18th 
(HOI; near-IR range). The near-IR model flux is scaled 
down by 15% compared to that in the optical, to al- 
low the simultaneous match to both wavelength regions. 
The model properties are the following: — 3.75 x 
lO^io, Tphot = 6260 K, i?phot = 12.9 x IQi^cm, Wphot = 
5530kms-\ pphot = 7.4 x lO-^^gcm-^, and n = 10. 
The model flux, which underestimates the observations 
by 10%, is renormalized to the observed value at 6IOOA. 
(This figure is available in color in the electronic version.) 

The overall fit quality is good. Discrepancies are no- 
ticeable for Ha, whose strength in absorption and emission 
is underestimated; we are currently investigating possible 
extra energy sources that could preserve the large spa- 
tial extent of Balmer/Paschen line formation, imposed by 
the morphology of Ha at late times. Contrary to previous 
epochs, the nitrogen abundance on this date has to be re- 
duced to maintain a good fit to the 8200A feature, whose 
strength is explained almost entirely by the Nai line (see 
below). We thus adopt N/He = lO^'^ (by number). The 
77OOA feature was at early-times associated with atmo- 
spheric extinction alone, although on the 14th and 19th 
of November, the absorption seemed to stretch further 
to the red, suggesting a contribution from the O i(3p-3s) 
triplet at 7775A. On the 5th of December 1999, follow- 
ing the reduced expansion velocity in the line formation 
region, one can unambiguously observe a P-Cygni line 
profile shape, blending with the atmospheric absorption 
only in the trough and blueward. As for earlier models, 
the reproduction of this line requires enhancing the oxy- 
gen abundance to 0/He = 10~^ (by number). However, 
this leads to an overestimated strength for the O i (3s-3p) 
triplet hue at 7984A; other O i hues at 8446A and 9264A 
overlap with Can, or Mgii and C i lines and are therefore 
not particularly useful as abundance constraints. We re- 
tain the same carbon abundance as in Sect. riTTl where 
the near-IR Ci feature at 1.07/xm required an enhance- 
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Fig. 8. Top: full (red) and continuum (blue) synthetic 
spectrum fits to observations (black) taken on the 5th of 
December 1999 (L02, day 37; this spectroscopic dataset 
is later used for distance determinations, in combination 
with the photometry on day 38 for which no spectra 
exist.) The model has the following parameters: = 
2 X 1O**L0, Tphot = 5920K, iiphot = 10.9 x lO^^cm, 
Wphot = 3400kms"\ pphot = 38.8 x 10~^''gcm~3, and 
n — 10. The model flux, which underestimates the ob- 
servations by 10%, is renormalized to the observed value 
at 7000A. Bottom: montage of rectified spectra computed 
by including bound-bound transitions of individual ioniza- 
tion species, ordered from the bottom by increasing atomic 
weight. For species labeled with an "a", flux departure 
from unity (continuum level) is magnified by a factor of 
four beyond 5680A. (This figure is available in color in the 
electronic version.) 



ment, with C/He = 1.7 x 10^"^. However, here, a clearly- 
defined C I feature is observed, stemming from the (3p-3s) 
multiplet covering the range 9060-9110A. We associate 
the 9200A feature with the Mgii (4p-4s) doublet at 9218- 
9244A, overlapping with the O i multiplet at 9260A (note 
that we adopt a solar magnesium abundance). Finally, the 
absorption strength of the Nai(3p-3s) doublet at 5890- 
5895A is well matched with a factor of four enhancement 
compared to solar (Prantzos et al. 1986). 



3.3. Discussion of Model Results 

The models presented in the previous section provide sat- 
isfactory fits to our sample of spectrophotometric obser- 
vations of SN1999em. Wc now investigate the overall ac- 
curacy of these model parameters, focusing on chemical 
abundances fSect. I!r3.1|l . outflow ionization/temperature 
and reddening (Sect. I3.3.2|l . and photospheric velocity 
(Sect. ESS. 

3.3.1. Chemical Abundances 

Using spectroscopic observations during the photospheric- 
phase evolution of type II SN, we anticipate inferring 
an outflow composition in conformity with CNO-cycle 
equilibrium values, typical of the envelope of the corre- 
sponding massive star progenitor, with at most traces 
of elements produced through explosive nucleosynthesis. 
Indeed, adopting a blue- or red-supergiant surface com- 
position (Prantzos et al. 1986) allows a good match to ob- 
served features of all species, for both their absolute and 
relative strengths. But let us now investigate how such line 
strengths are affected by a relative modulation of the out- 
flow H/He/CNO chemistry, as would result from a distinct 
evolutionary stage or chemical mixing efficiency (linked to, 
e.g., rotation and/or magnetic fields). 

We show in Fig. a comparison between the obser- 
vations of the 3rd of November 1999 (Sect. I3.2.3|l and 
two synthetic spectra corresponding to more (top) or 
less (bottom) evolved progenitor stars. For the former, 
we adopt Xh = 0.38, Xho = 0.60, Xc = 3.1 x lO""*, 
Xn = 1.4 X 10^^ and Xq = 2.4 x 10""^; for the lat- 
ter, we use Xu = 0.7, X^c = 0.28, Xc ^ 1.44 x lO^*, 
= 6.7 X 10-^, and Xq = 1.13 x 10"". 

The modulation in hydrogen abundance translates into 
a variation of the density of free-electrons in these fully 
ionized models: the higher the value of Xn, the higher 
the electron-scattering optical depth, either at the photo- 
sphere or at the model base. Here, we find that the pho- 
tosphere for the Xn — 0.7 model is located 5% further 
out than the Xn = 0.38 model, at a cooler temperature 
(10270K rather than 11 200 K), explaining the slightly 
"cooler" SED for the hydrogen-richer model. Beside this 
slight change in the slope of the SED, the strength of in- 
dividual features is noticeably altered. 

In the — 0.38 model, Balmer line strengths are sig- 
nificantly underestimated, and He 1 5875A overestimated. 
We also find that the increased nitrogen abundance leads 
to the overestimate of the strength of N ii features around 
5600A and 4600A. A comparison of Fig. M with Fig. El 
shows that the depleted hydrogen model is clearly of infe- 
rior quality compared to our best fit model. With a solar 
H/He ratio the changes relative to our best fit model are 
smaller, and it is more difficult to distinguish between the 
models. The sensitivity to changes in hydrogen, helium, 
and nitrogen composition is to be contrasted, for exam- 
ple, with Eastman & Kirshner (1989), who required an 



Luc Dessart and D. John Hillier: Type II supernovae spectroscopic modeling 



11 



unrealistically high hehum abundance to reproduce the 
Hel5875A hne profile in the early spectra of SN1987A. 

For carbon and oxygen, the situation is less clear: the 
outflow ionization is too low to lead to C ii or O il lines 
anywhere in the UV, optical, and near-infrared. At a stage 
where hydrogen recombines in the outflow, the ionization 
is suitable for the C i line to appear around 9000A, whose 
fitting requires an enhancement of about a factor of ten 
above that of the baseline model fSect. IX'2.8|I . Similarly, 
the same carbon abundance enhancement is needed to ex- 
plain the feature at 1.07/im, observed in the spectrum of 
the 19th of November (Fig.[7|). We find a similar situation 
for O I lines, which appear in the red part of the spectrum 
only at late times, and require a significant enhancement 
compared to CNO-cycle equilibrium values. Most notably, 
the Oi lines at TTOOA and 1.13/im require to have 0/He 
in the vicinity of 0.01 (by number). 

Overall, this investigation gives credence to our base- 
line model composition fSect. lXTjl at early times, but with 
C/0 abundances enhancing, as times progresses, above 
CNO-cycle equilibrium values; this suggests enhanced (ro- 
tational) mixing in the envelope of the progenitor star 
prior to collapse, or the appearance of the deeper and 
more chemically-evolved layers after just a few weeks past 
core-collapse. 

3.3.2. Outflow Ionization, Effective Temperature, and 
Reddening 

In the parameter space of photosphcric-phase type II SN 
(Paper I), the similar effects produced by temperature and 
reddening modulations on the synthetic SED motivate a 
joint discussion. 

In Fig. 1101 we present spectroscopic observations of 
SN1999em taken on the 5th of November 1999 fSect. lXTl 
and Fig. , but now overplotted with the SED from mod- 
els with decreasing photospheric temperatures: Tphot = 
9560 K (blue), Tphot = 9200 K (red), and Tphot = 8850 K 
(green) . Within the corresponding range of outflow ioniza- 
tion, the reduction of the UV flux is not so much caused by 
the decrease in Tphot, than by the opacity shift following 
the recombination from Fe III to Fe il in the SN outflow (as 
well as the opacity increase from numerous doubly-ionized 
metal species). In the optical, this shift has little bearing 
and the synthetic SED changes solely through the modest 
variation in Tphot- 

Note, however, that applying a reddening modulation 
to the model corresponding to the red line in Fig. 1111 
leads to a similar modulation: the relatively-enhanced UV- 
extinction (Cardelli et al. 1988) echoes the strong line- 
blanketing due to metals discussed above. The relatively- 
reduced optical-extinction echoes the weak SED sensitiv- 
ity to small Tphot variations. Overall, a change in E{B—V) 
of 0.05 or in Tphot of '^±500K have comparable impacts 
on the SED in the present parameter space. 

Early-epoch multi-wavelength observations are thus 
very important: metal line-blanketing is then only moder- 
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Fig. 9. Synthetic fits (red) to observations (black) of 
SN1999em on the 3rd of November 1999 (HOI; day 5) us- 
ing a model with the same parameters as those given in 
Sect. 13.231 fsee also Fig.Ol but having a more (top) or less 
(bottom) evolved chemical composition than the baseline 
model described in Sect. 13.11 (This figure is available in 
color in the electronic version.) 

ate and helps disentangle between temperature and red- 
dening effects; The UV range is also essential since this is 
where tests on the extinction magnitude are most sensi- 
tive. 

3.3.3. Photospheric velocity 

In Sect. 5 of Paper 2, we discussed how the velocity 
at maximum absorption in a sample of optical-range 
synthetic P-Cygni line profiles (i.e. Balmcr lines and 
Feii5169A) relates to the corresponding model photo- 
spheric velocity. Contrary to the usually held belief, such 
a measurement on optically-thick lines does not system- 
atically overestimate the model photospheric velocity. For 
any given ray, with impact parameter p, the location of 
maximum absorption for the flux-like quantity p ■ I{p) is 
exterior to the (continuum) photosphere, but for increas- 
ing p, this location shifts to smaller line-of-sight velocities, 
ultimately reaching line center. The total line flux is the 
integral of p ■ I{p) over all p, and thus reflects this range 
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SN1999em - 5lh of Nov. 1999 - E(B-V) = 0.1 
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Fig. 10. Illustration of the impact of a temperature 
change on the synthetic SED, implemented through a 20% 
increase (blue) or decrease (green) in luminosity (at con- 
stant base radius) compared to the model presented in 
Fig. El of Sect l3.2.4l (red) and observations (black). All 
other model parameters are held fixed. (This figure is 
available in color in the electronic version.) 



SN1999em - 5th of Nov. 1999 - E(B-V) = 0.1 




Fig. 11. Illustration of the effect of reddening modula- 
tions on the SED. We show the model and observations 
for the 5th of November 1999 (Fig. ^ day 7), reddening 
the synthetic SED with E{B-V) of 0.05 (blue), 0.1 (red), 
and 0.15 (green). Note how, in this parameter range, a 
0.05 change in E{B — V) is comparable to a ~±500K 
shift in photospheric temperature (Fig. I10|l . (This figure 
is available in color in the electronic version.) 



of velocities for the region of maximum-absorption. We 
find that although the velocity at maximum absorption 
for an optically-thick line always overestimates Wphot for 
the ray with impact parameter p~0, it can either overesti- 
mate or underestimate Wphot when the total line profile is 
considered. In practice, we identified few parameters that 
intervene in this correlation, the most important being the 
optical-depth of the line considered, which is related to the 
outflow ionization state, and the density gradient in the 



line formation region. These are physical entities that are 
difficult to constrain without detailed modeling. 

Second, P-Cygni line profiles have shapes that make 
the location of maximum absorption difficult to identify. 
This problem is severe at early times (say for the first 
week after explosion), when the expansion velocity in the 
continuum- and line-forming regions is much larger and 
the flux in the lines departs little from that in the contin- 
uum, making line features quite weak. This can be under- 
stood in the following way: suppose we take a model and 
increase its photospheric velocity: to first order, this will 
not change level populations, outflow ionization etc... and 
thus the total flux in all lines will be essentially the same, 
but the flux per-velocity-bin will be reduced, the more 
so for larger velocities. In practice, inspection of the first 
spectra for, e.g., SN1987A and SN1999em, show both very 
broad and very weak P-Cygni troughs, the region of maxi- 
mum absorption covering few thousand kms~^. In view of 
the argument discussed in the preceding paragraph, it is 
then difficult to infer a meaningful and accurate velocity 
for the photosphere. 

Third, the assumption of a unique radius for the pho- 
tosphere might not be appropriate on physical grounds. 
Multi-dimensional radiation-hydrodynamics simulations 
of core-collapse SN (Burrows et al. 1995; Janka & Miiller 
1996) predict explosions with a systematic departure from 
sphericity, both from the large-scale (following, e.g., the 
core and envelope rotation) and small-scale (following, 
e.g., Rayleigh- Taylor instabilities) viewpoints. Such effects 
are expected to make the photosphere quite "fuzzy" at 
later times (say one month after explosion), showing a 
range of spatial scales for different lines of sight and a 
significant dispersion in the expansion velocity. 

To illustrate the difficulty in assessing the photo- 
spheric velocity from P-Cygni line profiles, we show in 
Fig. El synthetic fits (red) to the observation (black) of 
SN1999em taken on the 9th of November 1999. The mod- 
els shown have identical parameters to the one presented 
in Sect. l3.231 departing only from its photospheric velocity 
Wphot.o = 7960 km s~^: if the model photospheric velocity 
is 10% higher (top), the Ha width is overestimated but 
H/3, Caii3800A, and Feii5159A are better fitted than in 
the model with a 10% slower photospheric velocity (bot- 
tom) . But all three models fit almost equally well the ob- 
servations so that there is a physically-meaningful uncer- 
tainty in photospheric velocity at the 10% level. While it is 
possible to measure precise velocities from the location of 
maximum absorption in some P Cygni profiles, the errors 
associated with such measurements significantly underes- 
timate the errors in the inferred photospheric velocity. 

Thus, the inference of the photospheric velocity is 
made difficult by the ambiguity of interpreting line width 
and velocity at maximum absorption in P-Cygni profiles, 
as well as, perhaps, the physical velocity range spanned 
by the line formation region of different species. 

For later use in distance determinations, we thus 
choose the model photospheric velocity obtained in 
Sect. 13. 21 and given in TableEl direct measurements on the 
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Fig. 12. Top: synthetic fits (red) to observations (black) 
of SN1999em on tlie 9tli of November 1999 (HOI; day 
11) using a model with the same parameters as those 
given in Sect. 13.2.51 (see also Fig. O except for a pho- 
tospheric velocity enhanced by 10% (wphot — 8760 km s~^, 
'^phot,o — 7960 km s^^). Bottom: same as above, but this 
time the photospheric velocity is reduced by 10% (wphot — 
7160 km s~^) compared to that of the reference model of 
Sect. 13.2.51 (This figure is available in color in the elec- 
tronic version.) 

observed spectra are more precise but the correspondence 
with t;phot is unclear. For all dates, the overall quality of 
our fits is good, which gives us confidence that our model 
photospheric velocity is a meaningful reflection of the ex- 
pansion velocity of the photosphere, say to within ±10% 
at most. 

4. The Distance to SN1999em 

Let us now determine the distance D to SN1999em — 
to do so we will use 5 distinct, but closely related, ap- 
proaches. For the flrst 2 approaches we use the "standard" 
form of the EPM (HOI, L02), in combination with either 
E96/H01's or Paper H's correction factors. In the third 
approach we use the models to derive the temperatures 
which are then combined with Papers IPs correction fac- 
tors. With the final EPM approach, we use both temper- 



atures and correction factors inferred from the best fit 
models. Finally, we use the SEAM approach of Baron et 
al. (2004). 

As discussed in Sects. I2I3I we selected eight epochs 
(October 30th, November 1st, 3rd, 5th, 9th, 14th, 19th, 
and December 5th) for which CMFGEN fitted well both 
continuum and line fiuxes. Relative to the discovery date 
HJD 2,451,480.94, these correspond to days 1, 3, 5, 7, 11, 
16, 21, and 38. No spectra were taken on day 38: for that 
date, given the very slow spectroscopic and photometric 
evolution at such late times, we combine the spectrum 
on day 37 with the photometry of day 38. Whatever the 
approach, we use the relevant model photospheric velocity. 
Note also that, in this section, all synthetic magnitudes 
are computed using the transmission functions and zero 
points of HOI. 

4.1. Inference with the EPM 

A key component of the EPM (E96, HOI, L02, Paper II) 
lies in the correspondence between correction factors and 
blackbody color temperatures Tc, defined from detailed 
models of type II SN spectra (E96; Paper II). In Paper 
II, we showed that our correction factors possess similar 
properties to those found in E96, albeit offset upwards by 
^20%. At a given Tc, correction factors show an intrinsic 
scatter of ^10-20%, roughly independent of Tc. However, 
if caused by uncertainties in this scatter increases dra- 
matically at low-Tc. Tabulated correction factors may also 
be inadequate to account for the blackbody assumption, 
if changes in magnitudes are due to variations in expan- 
sion velocity or metallicity, rather than changes, e.g., in 
temperature or the strength of flux dilution. Hence, we 
anticipate that using tailored models for each individual 
date will be more consistent and thereby improve the use 
of the EPM technique for distance determinations. 

In this section, we show a sequence of computations of 
distances based on distinct sets of quantities ^5, ^5 and 
T5. Here, S refers to one of the three bandpass combi- 
nations, {B,V}, {B,V,I}, and {V,I}, 9s to the quantity 
Rphot/D (not known individually), to correction fac- 
tors, and Ts to the blackbody color temperature, each of 
the last three quantities corresponding to a selected band- 
pass combination S. An error analysis of the EPM tech- 
nique and a description of the technique used to derive 
the distances are given in Appendix A & B. 

In Table El we show such a set computed with the 
"standard" form of the EPM (L02, HOI) and using 
HOl's correction factors. Restricting flrst the observational 
dataset to the first seven dates (prior to and excluding the 
5th of December 1999)^, we obtain increasing distances 

^ Throughout the paper we use observational data sets con- 
taining either 7 or 8 dates. This serves two purposes: First, Ha 
in the last data set is not well fitted, meaning that system- 
atic errors could, potentially, be large. Second, the effect of the 
longer time baseline, which helps to increase the accuracy of 
the explosion date, can be discerned. 
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for the three bandpass combinations, {B,V}, {B,V,I}, 
and {V,I}, from 8.6 Mpc, to 9.7 Mpc, and 11.7 Mpc, while 
the explosion date t^xp moves back in time in the same 
order, from -4.4 days, to -7.0 days, and -10.4 days. The 
disparity amongst bandpass combinations is large, espe- 
cially for the explosion dates, and unsatisfactory since the 
EPM should yield similar distances and explosion times 
whatever choice of bandpass combination, provided cor- 
rection factors across the range of color temperatures are 
appropriate. Adding the eighth date to the dataset leads to 
systematically smaller distances which are in better agree- 
ment but the explosion dates still differ substantially. The 
derived distance is approximately 8.7 Mpc which is 75% 
of the Cepheid-distance to NGC 1637 and SN1999em. 

We show in Table ^ the corresponding quantities ob- 
tained with Paper II's ^-prescription. These new correc- 
tion factors are indeed always higher than those shown in 
TableOlfor the same date, but the magnitude of the differ- 
ence varies from date to date. Hence, despite the preserved 
blackbody color temperatures between Table ful the non- 
uniform offset in correction factors between the two pre- 
scriptions induces a non-uniform offset in individual 6s- 
Resulting distances are systematically larger, using either 
7 or 8 dates, and the dispersion in the explosion times is 
reduced. 

We have investigated whether enforcing a fixed scal- 
ing a; > 1 on correction factors given by HOI would clarify 
these behaviors. As above, we find it introduces no notice- 
able modulation of the inferred color temperature (these 
remain constant at the one percent level). However, the 
individual values 9s obtained are then scaled by 1/x, ba- 
sically leaving ^5 x 9s constant. The resulting distances 
in that case are also scaled by x, but the explosion dates 
are now essentially unchanged. Using HOI or Paper II cor- 
rection factors thus modifies both the inferred time of ex- 
plosion and the distance, an effect arising from the non- 
uniform scaling in correction factors between HOI and 
Paper II; in other words, the differential variation of ^5 
with Ts affects distinct dates differentially. 

An important factor influencing the accuracy of 
the distance determinations is the incompatibility, at a 
given date, between color temperatures T5 obtained here 
(Table I3I4|I . and those previously computed on corre- 
sponding CMFGEN model SEDs (Table H in Sect.l^. 
For bandpass combinations {B,V,I} and {V,I}, model 
values of Ts are systematically larger than those deduced 
from the minimization technique performed on observa- 
tions, while for {B, V}, they are smaller at high and larger 
at low Tbv- This color- mismatch is likely to stem from 
the flux offset between synthetic and observed spectra, 
leading to an over- or under-estimate of the correspond- 
ing bandpass magnitudes on which the computation of 
color-temperatures is based. Figure |3| shows such a sit- 
uation where the observed / band flux is substantially 
underestimated by our model - note that the flux cali- 
bration may also be at fault, especially at longer optical 
wavelengths where atmospheric extinction operates. To be 
more quantitative, we provide an estimate of the effect of 



flux errors on the derived angular diameter in Appendix 
A (Table A.l). 
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Fig. 13. Linear fit to the BV data shown in Table for 
8 data points (solid blue line) and 7 data points (bro- 
ken green line). The reciprocal of the slope is the dis- 
tance to the SN. The time is measured with respect to JD 
2451480.9. The vertical bars indicate ±10%. (This figure 
is available in color in the electronic version.) 

We now compute the EPM-distance using the model 
color-temperatures. First, we use the ^-prescription from 
Paper II and show resulting EPM quantities in Table |31 
Distances derived using the first 7 dates are systematically 
higher than the Cepheid distance, but have large errors. 
The explosion dates, within errors, are consistent. With 8 
dates, all distances agree within 1 sigma of the Cepheid 
distance of 11.7Mpc (L03) and the SEAM distance of 
12.4 Mpc (B04, see next section). Because the observa- 
tional and model color temperatures show filter set and 
epoch dependent differences, the changes in the derived 
distances between those shown in Table 0] and Table |5l are 
not systematic or uniform. 

Finally, combining both color-temperatures and cor- 
rection factors from the models of Sect. 13.21 we find 
even better agreement in the distance predicted between 
the different band passes, of 11.7Mpc, 11.6Mpc, and 
11.5 Mpc, in the usual sequence (TableE)). This is expected 
since the same models are used — indeed the differences, 
typically less than 2%, must simply refiect errors in inter- 
polations etc arising in the application of the EPM tech- 
nique. The corresponding explosion dates are -6.2 days, 
for all three data sets. In Fig. ^1 we show a linear fit of 
the form 9/vp\iot = {t — texp)/D to the BV data shown in 
Table El 

4.2. Inference based on CMFGEN synthetic SEDs 

The Spectral-fitting Expanding Atmosphere Method 
(SEAM; B04; Mitchefl et al. 2002) is an attractive alter- 
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Table 3. Table of EPM quantities 9s (10* kmMpc"-^), £_s, and Tg (K), for the eight epochs covered by our analysis 
(i.e. from the 30th of October till the 5th of December 1999; see Sect. Also given are the resulting distance and 
explosion dates inferred. Days and dates are relative to JD 2,451,480.94 



EPM combined with Minimization Technique and ^ prescription from HOI 



Day 


Angular size 


Correction Factor 


Temperature 




(10** km Mpc~ 












(K) 






Obv 


Obvi 


Ovi 


£,BV 


£,BVI 




Tbv 


Tbvi 


Tvi 


1.0 


6.28 


7.78 


8.84 


0.405 


0.425 


0.440 


16860 


13513 


11543 


3.0 


6.73 


8.31 


9.37 


0.401 


0.423 


0.440 


16476 


13241 


11383 


5.0 


8.16 


9.14 


9.82 


0.382 


0.416 


0.438 


14394 


12312 


10998 


7.0 


8.46 


9.89 


10.87 


0.377 


0.413 


0.436 


13801 


11479 


10070 


11.0 


11.10 


11.48 


11.85 


0.367 


0.412 


0.436 


11271 


10214 


9477 


16.0 


15.77 


13.89 


12.23 


0.461 


0.434 


0.437 


7747 


8532 


9365 


21.0 


16.45 


15.47 


13.89 


0.600 


0.478 


0.445 


6514 


7411 


8404 


38.0 


17.51 


17.58 


16.98 


1.202 


0.652 


0.482 


4688 


5777 


7043 






BV Set 






BVI Set 






VI Set 




Using the first 7 dates only 


D 




8.6 ±0.8 






9.7 ± 1.0 






11.7± 1.5 




icxp 




-4.4 ±0.9 






-7.0 ± 1.4 






-10.4 ±2.1 




Using 8 dates 


D 




8.0 ±0.6 






8.8 ±0.7 






10.1 ±0.9 




icxp 




-3.9 ±0.7 






-5.9 ± 1.0 






-8.4 ± 1.4 





Table 4. Same as Table El this time using Paper IPs ^-prescription 



EPM combined with Minimization Technique and ^ prescription from Paper II 



Day 


Angular size 


Correction Factor 


Temperature 




(10** kmMpc" 












(K) 






Obv 


Obv I 


Ovi 




£,BVI 


ivi 


Tbv 


Tbvi 


Tvi 


1.0 


5.82 


6.58 


7.41 


0.436 


0.504 


0.526 


16876 


13481 


11543 


3.0 


6.20 


6.95 


7.86 


0.439 


0.505 


0.525 


16364 


13241 


11351 


5.0 


6.88 


7.48 


8.16 


0.453 


0.508 


0.525 


14378 


12328 


11031 


7.0 


6.95 


7.93 


8.99 


0.460 


0.514 


0.526 


13769 


11495 


10086 


11.0 


8.09 


8.99 


9.74 


0.503 


0.530 


0.530 


11271 


10166 


9493 


16.0 


10.57 


10.50 


10.12 


0.688 


0.572 


0.532 


7747 


8548 


9317 


21.0 


11.63 


11.70 


11.25 


0.854 


0.632 


0.547 


6498 


7411 


8420 


38.0 


14.87 


13.96 


13.51 


1.414 


0.820 


0.606 


4688 


5777 


7043 






BV Set 






BVI Set 






VI Set 




Using the first 7 dates only 


D 




13.5 ± 1.5 






12.5 ± 1.6 






14.6 ± 1.9 




^cxp 




-7.2 ± 1.4 






-8.4 ± 1.7 






-10.1 ± 2.2 




Using 8 dates 


D 




11.7± 1.0 






11.9 ± 1.0 






12.6 ± 1.2 




^cxp 




-5.7 ± 1.0 






-6.9 ± 1.2 






-8.8 ± 1.5 





native to the EPM: each observed SED is directly fitted 
with a synthetic SED, produced by a detailed model at- 
mosphere calculation for the corresponding date, avoiding 
necessary corrections for approximating observed colors 
with those of a single-temperature blackbody (EPM; see 
previous section). 

Nonetheless, the SEAM follows a similar procedure to 
the EPM. Indeed we would expect to derive very simi- 
lar answers to that found using the EPM technique with 
detailed SN fitting. We first select a number of epochs 



for which photometric and spectroscopic observations are 
available and accurate, and for which we possess high qual- 
ity fits to the observed SED. We then select a range of 
explosion times, and compute the photospheric radius ex- 
pected for each date of observation, given by -Rphot(i) = 
'fphot(i) (i — ^exp), where t is the observation date, icxp the 
explosion date, and Wphot(i) represents the corresponding 
photospheric velocity (the initial radius of each considered 
mass shell is neglected). As in the EPM, the assumption 
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Table 5. Same as Table 13 this time using blackbody color temperatures from Sect.O and Paper IPs ^-prescription 



EPM combined with Minimization Technique, 
blackbody color temperatures from Sect. El 
and ^-prescription from Paper II 



Day 


Angular size 


Correction Factor 


Temperature 




(10** kmMpc' 












(K) 






Obv 


Obvi 


Ovi 


^BV 


£,BVI 




Tbv 


Tbvi 


Tvi 


1.0 


6.13 


5.90 


5.67 


0.440 


0.503 


0.543 


16024 


14965 


14388 


3.0 


6.35 


6.05 


5.97 


0.441 


0.503 


0.542 


15959 


15030 


14228 


5.0 


6.95 


6.95 


7.03 


0.455 


0.505 


0.530 


14228 


13202 


12529 


7.0 


7.33 


7.63 


7.86 


0.470 


0.510 


0.525 


12945 


11951 


11278 


11.0 


8.46 


8.91 


9.37 


0.519 


0.529 


0.527 


10701 


10252 


9835 


16.0 


10.12 


10.12 


9.82 


0.650 


0.562 


0.529 


8168 


8841 


9547 


21.0 


10.72 


11.03 


10.80 


0.760 


0.605 


0.540 


7110 


7847 


8745 


38.0 


12.08 


12.76 


13.21 


1.077 


0.734 


0.594 


5539 


6340 


7238 






BV Set 






BVI Set 






VI Set 




Using the first 7 dates only 


D 




14.2 ± 1.6 






13.2 ± 1.4 






12.9 ± 1.3 








-8.2 ± 1.6 






-7.2 ± 1.4 






-6.9 ± 1.3 




Using 8 dates 


D 




12.8 ± 1.1 






12.0 ± 1.0 






11.7 ±0.9 




^cxp 




-8.0 ± 1.2 






-6.2 ± 1. 






-5.9 ± 1.0 





Table 6. Same as Tabled this time using blackbody color temperatures and ^-values from Sect. 13 



EPM combined with blackbody color temperatures 
and 5- values from Sect.|3 



Date 


Angular size 


Correction Factor 


Temperature 




(10* kmMpc" 












(K) 






6bv 


Obvi 


Ovi 


^BV 


£,BVI 


ivi 


Tbv 


Tbvi 


Tvi 


1.0 


5.82 


5.90 


5.97 


0.462 


0.500 


0.519 


16024 


14965 


14388 


3.0 


6.50 


6.58 


6.65 


0.430 


0.462 


0.488 


15959 


15030 


14228 


5.0 


7.11 


7.18 


7.26 


0.443 


0.488 


0.513 


14228 


13202 


12529 


7.0 


7.78 


7.86 


7.86 


0.443 


0.494 


0.526 


12945 


11951 


11278 


11.0 


9.37 


9.44 


9.37 


0.469 


0.500 


0.526 


10701 


10252 


9835 


16.0 


10.20 


10.20 


10.27 


0.646 


0.558 


0.507 


8168 


8841 


9547 


21.0 


11.85 


12.01 


12.16 


0.685 


0.558 


0.481 


7110 


7847 


8745 


38.0 


12.16 


12.38 


12.83 


1.072 


0.755 


0.608 


5539 


6340 


7238 






BV Set 






BVI Set 






VI Set 




Using the first 7 dates only 


D 




12.4 ± 1.2 






12.4 ± 1.3 






12.4 ± 1.3 




^cxp 




-6.8 ± 1.3 






-6.9 ± 1. 






-7.0 ± 1.3 




Using 8 dates 


D 




11.7 ± 1.0 






11.6 ± 1.0 






11.5 ±0.9 




icxp 




-6.2 ± 1.0 






-6.2 ± 1.0 






-6.2 ± 1.0 





of homologous expansion, warranted past day after 
core-bounce, is essential. 

In our approach, we do not evolve an original hydrody- 
namical input for the exploding envelope; we use instead 
an analytical description of the density distribution, and 
adjust the luminosity to match typical values expected 
from explosion models and former spectroscopic studies 
(for SN1999em, we use the Cepheid-distance to ensure 
our adopted luminosities are compatible with observed 



fluxes and magnitudes). For a chosen iexp corresponds an 
^phot(Oi usually different from that of the fitting model 
for that date; we thus scale the model luminosity in pro- 
portion with the square of the change in i?phot(Oj ^-iid 
update synthetic magnitudes accordingly. Then, for each 
choice of the explosion date texp and for each observation 
epoch we compute the distance modulus nsitexp,t) for 
all bandpass sets S — {B, V}, {B, V, /}, and {V, /}, given 
by Ms = (a*x)s, where ^ix = mx - Mx - Ax; X repre- 
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sents one of the band passes in the set S; mx, Mx, and 
Ax are the apparent, absolute, and extinction magnitudes 
in bandpass X. Thus, for each chosen explosion date, we 
compute, for each set 5, a distance modulus ^s(ioxpj t) av- 
eraged over all observation epochs. The distance modulus 
selected at the end corresponds to the ^cxp that minimizes 
the scatter of ns{texp,t) over all epochs. 

We show results for this sample of observations in 
Table [3 The results are grouped by triads. Using the 
adopted reddening E{B — V) = 0.1 and only seven 
(eight) observations, we obtain a distance of ^12.15 Mpc 
(11.5 Mpc), in good agreement with the SEAM-distance of 
B04 (12.4 ±2.5 Mpc) or the Cepheid distance of Leonard 
et al. (2003; 11.7 ± 1 Mpc). Limiting the sample to the first 
four epochs only leads to a significant increase in the dis- 
tance, to ~15.2Mpc: the much reduced time-span makes 
this prediction unreliable. Note also that the use of the 
.Z— band has only a modest impact on the resulting dis- 
tance. 

If we assume that our adopted models, which fit well 
the observations with the choice of E{B — V) — 0.1, 
would still fit the observations adequately with a mod- 
ulation of the reddening by ±0.05 about this nominal 
value (see the discussion in Sect. 13.3"^ . then such a red- 
dening uncertainty introduces a distance uncertainty of 
« ±0.8 Mpc (equivalent to a ±0.15 range in distance 
modulus). However, this represents an upper limit since 
new models should be computed that fit better the obser- 
vations with the newly adopted reddenings. When this is 
done, the change in derived temperatures and radii par- 
tially compensate for the change in reddening (Appendix 
A). 

We find that a ± 10% error in Wphot has a negligi- 
ble impact on the inferred time of explosion (constrained 
by the relative evolution of magnitudes between differ- 
ent epochs) but affects noticeably the inferred distance: 
at a given epoch, smaller (larger) photospheric veloci- 
ties mean smaller (larger) photospheric radii and thus 
smaller (larger) distances to match the observed magni- 
tudes. Thus, by itself, the error on fphot fSect. l^3.3() intro- 
duces an uncertainty of ± 1 Mpc in the inferred distance. 
Finally, we have added an eighth observation to our sam- 
ple used so far (that of the 5th of December 1999). The 
average distance for the three bands used is then reduced 
by ~0. 5-1 Mpc compared to that of the reference case. 

5. Conclusion 

In this paper, we have used the model atmosphere code 
CMFGEN (Paper I and II) to perform a quantitative spec- 
troscopic analysis of the photospheric-phase evolution of 
SN1999em, a type II supernovae. 

We obtain high-quality fits to fiux-calibrated data 
for eight epochs (from the 30th of October until the 
5th of December 1999), modeling the evolution of the 
optically-thick outfiow from fuU-ionization to recombina- 
tion. Indeed, we find that the spectroscopic changes are 
primarily due to the temperature decrease following ex- 



Table 7. SEAM results for the explosion date and dis- 
tance to SN1999em, under a variety of choices that high- 
light uncertainties: bandpass set, number of observation 
dates, reddening, and photospheric velocity values. Time 
is relative to JD 2,451,480.54. 
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pansion and radiation-leakage, associated with a decrease 
in the outflow ionization. The critical shift from Fciii to 
Fell at Tphot ^ 8000 K, concomitant with the onset of 
hydrogen recombination, causes a dramatic increase in 
opacity, first in the UV and later in the optical, that 
alters considerably the overall shape of the spectral en- 
ergy distribution. For fully or partially ionized models, we 
were able to reproduce, with good accuracy, the shape 
of both the spectral energy distribution and line profiles. 
The thus-constrained chemical composition is compatible 
with CNO-cycle equilibrium abundance values for the first 
month. At early times, the carbon and oxygen abundances 
are poorly constrained due to the lack of lines of Cii and 
On in the optical and near-IR. At later times, with fur- 
ther cooling of the outflow, C i and O i lines are seen in the 
7000-12 OOOA regions, requiring signiflcant enhancements 
over the CNO-cycle equilibrium values; it is not clear if 
this corresponds to a genuine increase in the correspond- 
ing abundances as the photosphere recedes into deeper 
layers, or if this abundance pattern was identical at all 
previous times. 

Significantly, there was no observational incentive 
to vary the metal composition, whatever the epoch, 
so we adopted a unique mixture corresponding to so- 
lar. In other words, over a 38-day period, we see 
no spectroscopic evidence of products from explosive- 
nucleosynthesis. Additionally, we find no evidence for a 
flatter density distribution at later times: for all epochs, 
we adopt a density exponent n — 10—12 (Paper I- 
II) which permits satisfactory fits to line profile shapes. 
Finally, as for other groups (e.g. Mitchell et al. 2002), 
we have difficulty in reproducing Balmer (and Paschen, 
when available) line profiles, both in strength and width. 
This discrepancy might stem from a missing outflow en- 
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ergy source, perhaps, e.g., in the form of energy deposition 
from radioactive decay of ^^Ni (Mitchell et al. 2002), or 
result from the neglect of time-dependence effects (see, 
e.g., Utrobin & Chugai 2005) in standard radiative trans- 
fer modeling of photospheric phase type II SN. We are 
currently studying this issue. 

We also estimate the distance to SN1999em using a 
variety of approaches and compare with the reference 
value of Leonard et al. (2003) using NGC 1637 Cepheids. 
Although we reproduce the under-estimated EPM dis- 
tance of Leonard et al. (2002) or Hamuy et al. (2001), there 
are non-trivial differences between our approach and as- 
sumptions, and these two works. First, our adopted photo- 
spheric velocities are higher or lower than theirs, depend- 
ing on the epoch: systematically higher photospheric ve- 
locities lead to faster increasing photospheres, which under 
the same conditions of magnitude etc., require a higher 
distance - the difficulty here is that the difference in veloc- 
ity is positive or negative at different dates. Second, the 
latest epoch is limited to day 21 (7th observation in our 
sample) or day 38 (8th observation) after discovery - add 
^^6.2 days to time of core-collapse - in our approach rather 
than 70days: there is a sensible incentive in lengthen- 
ing the time baseline to reduce the errors (see Appendix), 
but as argued in Paper II, the corrupting effect of lines 
in all optical band passes as well as the likely energy con- 
tribution from radioactive decay make such late time ob- 
servations less adequate for the EPM. Third, to extract 
the distance and time of explosion from our collection of 
{t, 0/v} data points, we use time as the independent vari- 
able, rather than 0/v: this seems a better approach since 
time is accurately known. Moreover, in our least-square 
fitting, we adopt constant fractional errors in 0/v rather 
than constant weights. Wc find that this approach yields 
distances that are systematically and at least 10% greater 
that those obtained by assuming d/v as the independent 
variable and constant weights. 

Applying the correction-factor recipe of Paper II leads 
to a higher distance, in good agreement with the Cepheid 
estimate. It thus seems that correction factors are indeed 
a major uncertainty in the EPM, the systematically lower 
values used previously playing a central role in the dis- 
crepancy (see, e.g.. Baron et al. 2004). Overall, however, 
approaches based on tabulated correction factors lead to 
a significant disparity in distances obtained from different 
bandpass sets, in conflict with the near-constancy, at early 
times, of the angular size of the photosphere with optical- 
wavelength. The problem likely stems from the intrinsic 
scatter of correction factors at a given color-temperature, 
obtained from observations alone or even adopted from 
tailored models. 

When both color-temperatures and correction factors 
are adopted from individual best-fit models to each ob- 
servation, the EPM technique delivers more consistent 
distance estimates between band passes. Slight differ- 
ences persist from small but systematic discrepancies be- 
tween the model and observed spectral distributions. But 
with such an approach, we obtain an EPM-distance of 



11.5±lMpc, in agreement with the Cepheid estimate. 
Hence, while there is nothing fundamentally wrong with 
the EPM, the need for detailed model atmosphere calcu- 
lations defeats the original purpose of the method. 

We obtain another distance estimate of 12.2±2Mpc 
using synthetic spectra and magnitudes, minimizing the 
scatter in individual distance moduli for a range of 
adopted explosion dates. This approach is similar to the 
approach of Baron et al. (2004) and delivers a similar dis- 
tance. The inferred explosion date is the 23rd of October 
1999 with an uncertainty of 1 day — this is approximately 
4 days earlier than the estimate of Hamuy et al. (2001) but 
similar to that of Leonard et al. (2002) and Baron et al. 
(2004). 

This study shows that there are a variety of reliable 
ways to determine distances with early-time photospheric- 
phase type II SN, with good prospects for cosmology. In 
forthcoming studies, we will reproduce the present analy- 
sis with a larger sample of objects, covering a wider range 
of distances or redshifts, and thereby attempt to provide 
an alternative, to using type la SN, for the cosmological 
distance scale. 
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Appendix A: Errors 

As an accurate distance is the goal of the EPM technique 
it is worth discussing the error sources, and how they effect 
the determination of the angular diameter, and hence the 
distance, in more detail. For simplicity we consider the use 
of two filters only. Let Fi be the observed flux in bandpass 
i, and the reddening corrected flux. Thus 



F^ = F, exp {QmiR,E{i ~ j)) 



(A.l) 



where Ri — Ai/E{i — j) and E{i — j) is the color ex- 
cess (e.g., E{B — V)). The blackbody color temperature 
is determined from 



#exp(0.921^;(.-,))=f^^ ^^P^"^)-' 
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exp(ui) - 1 



(A.2) 



where ui = hvi/KT. Under the assumption that the errors 
are small, we can determine them using standard tech- 
niques (Taylor 1997). Differentiating Fan. IA.2l vields 
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where 
G, = 



Table A.l. Errors in 9 due to flux and reddening errors. 



1 - exp(~Ui) 
The expression for the angular diameter, 9 is 



(A.5) 



(A.6) 



To allow for the error in ^ we split it into two parts, 
a random error at fixed temperature, ct^t, and an error 
attributed to the error in T which can be attributed to an 
error in Fi, Fj, and E{i — j). Thus, for example, 

_ d£, dT 
dFi ^ dfm 

and where d£,/dT is determined from a set of models. 
If we assume the errors are independent we obtain 
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The value of u^t is difficult to determine. An estimate 
of its value can easily be determined from a grid of models 
at fixed effective temperature encompassing a range of SN 
type models (e.g., density exponent). However, there is 
no guarantee that this grid is representative of true SN. 
A further source of error is the fact that models do not 
exactly match the observations, and this will contribute an 
error greater than those introduced by the dispersion in 
SN properties, This error could, if known, be incorporated 
into tr^T- 

In Table 10 we list G,/(G, - G^), d\n9/d\nF, 
dln9/dE{B — V) for three different passband combina- 
tions. We ignore the coupling with ^, and assume that the 
fractional errors in the two pass bands are identical. As 
readily apparent, the fractional error in 9 due to an error 
in flux increases with SN temperature — this is simply 
due to the increasing insensitivity of the blackbody slope 
to temperature as we approach the Rayleigh-Jeans regime. 
We also note that a 1% error in flux is generally more im- 
portant than a 0.01 error in E(B-V). This somewhat sur- 
prising result occurs since the effect of the reddening on 
the temperature estimate, and the direct effect of redden- 
ing on the angular diameter (because of a change in flux) 
partially compensate. This has been previously noted by 
Leonard et al. (2002), Eastman et al. (1996), and Schmidt 
et al. (1992). Not surprisingly, we also see that the BI 
passband combination, with its wider wavelength cover- 
age, gives potentially the smallest errors. Finally we note 
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that in the BV band, the angular diameter is very sensi- 
tive to flux errors. A random error of 3% in each passband 
leads to a 10 to 14% error in 9. 

Errors associated with ^ are specific to the EPM tech- 
nique. However the other error sources will also apply, in 
a complicated fashion, to the SEAM technique. The pri- 
mary distinction is that the SEAM technique uses global 
fitting, whereas the EPM is sensitive to local errors in 2 
or 3 pass bands. These errors arise, for example, from er- 
rors in the abundance, the neglect of some species, from 
observational errors, and from inadequacies in the model. 

Appendix B: Distance determination 

The relation between the distance and 9 is 

9 _ (t — <oxp) 



D 



(B.l) 



Several different approaches can be used to estimate D 
and tcxp- Because t is known accurately, we must use t as 
the independent variable. In standard form the equation 
is 



y ^ a 



bt 



(B.2) 



If we use least squares, errors for a and b can be readily 
determined from the accuracy of the fit. Thus we have 
D = 1/6 and 



<7d/D = ab/b 
Similarly, t^xp 



,2 2 



a/b and 
2t 



(B.3) 



(B.4) 



where aab is the covariance of a and b. 

The errors on 9/v are difficult to determine accurately. 
However, in general we expect the errors in 9/v to scale 
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with 9 /v. For simplicity, we therefore perform weighted 
least squares fitting assuming constant fractional errors. 
Closer examination of the error analysis discussed earlier 
suggests that the fractional errors in 0/v could be some- 
what larger at earlier epochs, but given the additional 
influence of systematic errors, we have made no effort to 
take this into account. 

We also used robust fitting techniques (see Press et al. 
1986) with weights, but the answers were generally similar 
(i.e., within 1 sigma) to those obtained with least squares 
fitting. As a test, we also tried fitting the data assuming 
constant errors. However, this gives too much weight to 
the data of day 38, and leads to unrealistically small error 
bars. 

The errors on D and <cxp listed in the tables are from 
the fit only, and thus are a MINIMUM estimate of the er- 
ror. No systematic errors, which are difficult to determine, 
have been taken into account. Systematic errors will arise 
from errors in reddening, and possibly, errors in the model. 
It is also likely, for example, that errors in adjacent 6/v 
estimates arc correlated, since the spectra are similar, and 
hence model errors are likely to be similar. 
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